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1 Inline Sage, code blocks 

This is an example 2 + 2 = 4. If you raise the current year mod 100 (which 
equals 15) to the power of the current day (17), you get 98526125335693359375. 
Also, 2015 modulo 42 is 41. 

Code block which uses a variable s to store the solutions: 

1+1 

var ( ' a,b, c ' ) 

eqn = [a+b*c==l, b-a*c==0, a+b==5] 
s = solve(eqn, a,b,c) 

Solutions of eqn = [be + a = 1, — ac + b = 0, a + b = 5]: 

25i\/79 + 25 , hi V79 + 5 1 /— 1 

a = 1= , b = — , c = — z v 79 + — 

6i\/79 - 34 iV79 + U 10 10 

25^^79 - 25 , 5iV79-$ 1 1 

a = 7 = , b = — — = , c = -IV79+ — 

6iV79 + M i \/79 - 11 10 10 

Now we evaluate the following block: 

E = EllipticCurve("37a") 

You can't do assignment inside \sage macros, since Sage doesn't know how to 
typeset the output of such a thing. So you have to use a code block. The elliptic 
curve E given by y 2 + y = x 3 — x has discriminant 37. 

You can do anything in a code block that you can do in Sage and/or Python. 
Here we save an elliptic curve into a file. 

try : 

E = load('E2') 
except IDError: 

E = EllipticCurveC [1,2,3,4,5] ) 

E.anlist(100000) 

E.save('E2') 
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The 9999th Fourier coefficient of y 2 + xy + 3y = x 3 + 2x 2 + Ax + 5 is -27. 
The following code block doesn't appear in the typeset file. . . but we can 
refer to whatever we did in that code block: e = 7. 

var ( ' x ' ) 

f(x) = log(sin(x)/x) 

The Taylor Series of / begins: x — 467 1 775 x 10 — 37 g 00 x 8 — ^"ilb xA ~ \ x ' ''■ 

2 Plotting 

Here's a very large plot of the elliptic curve E. 
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You can use variables to hold plot objects and do stuff with them. 

p = plot(f , x, -5, 5) 
Here's a small plot of / from —5 to 5, which I've centered: 



On second thought, use a size of 3/4 the \textwidth and don't use axes: 



Remember, you're using Sage, and can therefore call upon any of the software 
packages Sage is built out of. 

f = maxima ( ' sin(x) ~2*exp(x) ' ) 
g = f . integrate ( 'x' ) 

Plot g(x), but don't typeset it. 

You can specify a file format and options for includegraphics. The default 
is for EPS and PDF files, which are the best choice in almost all situations. 
(Although see the section on 3D plotting.) 
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If you use regular latex to make a DVI file, you'll see a box, because DVI 
files can't include PNG files. If you use pdf latex that will work. See the 
documentation for details. 

When using \sageplot, you can pass in just about anything that Sage can 
call .saveO on to produce a graphics file: 
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To fiddle with aspect ratio, first save the plot object: 



p = plot(x, 0, 1) + circle((0,0) , 1) 
p . set_aspect_ratio(l) 



Now plot it and see the circular circle and nice 45 degree angle: 
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Indentation and so on works fine. 

s =7 

s2 = 2~s 

P.<x> = GF(2) [] 

M = matrix (parent (x) ,s2) 

for i in range (s2): 

p = (l+x)~i 

pc = p.coeffsO 

a = pc. count (1) 

for j in range (a) : 

idx = pc. index (1) 

M[i,idx+j] = pc. pop (idx) 

matrixprogram = matrix_plot (M, cmap=' Greys ' ) 

And here's the picture: 



0 20 40 60 80 100 120 
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Reset x in Sage so that it's not a generator for the polynomial ring: x 



2.1 Plotting (combinatorial) graphs with TikZ 



Sage now includes some nice support for plotting graphs using TikZ Here, we 
mean things with vertices and edges, not graphs of a function of one or two 
variables. 

The graphics in this section depends on the tkz-berge package, which is 
generally only available in newer TjrpC distributions (for example, TgXLive 2011 
and newer). That package depends in turn on TikZ 2.0, which is also only 
available in newer TgX distributions. Installing both of those is in some cases 
nontrivial, so this section is disabled by default. 

If you have TikZ and tkz-berge and friends, remove the comment environ- 
ments below. 



2.2 3D plotting 

3D plotting right now (Sage version 4.3.4) is problematic because there's no 
convenient way to produce vector graphics. We can make PNGs, though, so 
if you pass sageplot a graphics object that cannot be saved to EPS or PDF 
format, we will automatically save to a PNG file, which can be used when 
typesetting a PDF file, but not when creating a DVI file. However, you can 
specify the "imagemagick" option, which will use the Imagemagick convert 
utility to make EPS files. See the documentation for details. 

Here's a 3D plot whose format we do not specify; it will automatically get 
saved as a PNG file and won't work when using latex to make a DVI file. 
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3 Pausing SageTgX 



Sometimes you want to "pause" for a bit while writing your document if you 
have embedded a long calculation or just want to concentrate on the M^]X and 
ignore any Sage stuff. You can use the \sagetexpause and \sagetexunpause 
macros to do that. 

A calculation: (SageTJrX is paused) and a code environment that simulates 
a time-consuming calculation. While paused, this will get skipped over. 

import time 
time . sleep(15) 



Graphics are also skipped: 



SageTfrX is paused; no graphic 



4 Make Sage write your BT^X for you 

With SageT^X, you can not only have Sage do your math for you, it can write 
parts of your ffl^X document for you! For example, I hate writing tabular envi- 
ronments; there's too many fiddly little bits of punctuation and whatnot. . . and 
what if you want to add a column? It's a pain — or rather, it was a pain. Just 
write a Sage/Python function that outputs a string of M^X code, and use 
\sagestr. Here's how to make Pascal's triangle. 

def pascals_triangle (n) : 

# start of the table 

s = [r"\begin{tabular}{cc| " + "r" * (n+1) + "}"] 
s.appendCr" & & $k$: & \\") 

# second row, with k values: 
s.appendCr" & ") 

for k in [0. .n] : 

s. append ("& {0} ".format(k)) 
s . append (r"\\") 

# the n = 0 row: 

s.append(r"\hline" + "\n" + r"$n$: & 0 & 1 & \\") 

# now the rest of the rows 
for r in [1 . .n] : 

s. append (" & {0} ".format(r)) 
for k in [0. .r] : 

s. append ("& {0} " .format (binomial (r, k) ) ) 
s. append (r"\\") 

# add the last line and return 
s . append(r"\end{tabular}") 
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return ' ' . join(s) 

# how big should the table be? 
n = 8 

Okay, now here's the table. To change the size, edit n above. If you have 
several tables, you can use this to get them all the same size, while changing 
only one thing. 
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5 Include doctest-like examples in your docu- 
ment 

Here are some examples of using the sageexample environment: 
sage: 1+1 

2 

sage: factor (x~2 + 2*x + 1) 

(x + 1) 2 

If you want to see the plain-text output as well as the typeset output, renew 
the sageexampleincludetextoutput command to True: 

\renewcommand{\sageexampleincludetextoutput}{True> 

This can be useful to check that the two outputs are consistent. 

When this environment is near the bottom of the page, it may look like the 
page number is the output of a command, when in fact the real output is on the 
next page. If the output of a command below looks like 11, don't worry, that's 
just the page number. 

sage: 1+1 
2 
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sage: factor (x~2 + 2*x + 1) 
(x + 1)~2 

(x + lf 

Multiline statements are supported, as are triple-quoted strings delimited by 
single quotes: 

sage : def f (a) : 

... '''This function is really quite nice, 

although perhaps not very useful . ' ' ' 
print "f called with a = ", a 

y = integrate (SR(cyclotomic_polynomial (10)) + a, x) 
return y + 1 

sage: f (x) 

- x 5 - - x 4 + - x 3 + X + 1 
5 4 3 

Note that the "/ called with. . . " stuff doesn't get typeset, since when running 
Sage on example . sagetex . sage, that gets printed to the terminal. 

When typesetting your document, the validity of the outputs is not checked. 
In fact, the provided outputs are completely ignored: 

sage: is_prime(57) 
toothpaste 

False 

However, typesetting your document produces a file named example_doctest . sage 
containing all the doctest-like examples, and you can have Sage check them for 
you with: 

$ sage -t example_doctest . sage 

You should get one doctest failure from the "toothpaste" line above. 
Please look into this file for the original line numbers. 

Beware that sage -t does not handle well file names with special charac- 
ters in them, particularly dashes, dots, and spaces — this ultimately comes from 
the way Python interprets import statements. Also, running doctests on files 
outside the main Sage library does not always work, so contact sage-support 
if you run into troubles. 

Some more examples. This environment is implemented a little bit differently 
than the other environments, so it's good to make sure that definitions are 
preserved across multiple uses. This will correctly define a, but not print its 
output because the statement is made up of a sequence of expressions. 
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sage: 1; 2; a=4; 3; a 

After that, Sage should remember that a = 4 and be able to use that in future 
sageexample blocks: 

sage: f(a) 

1 5 1 4 1 3 1 2 r 

- x b x i + -x fi x 2 + 5 x + 1 

5 4 3 2 

6 Plotting functions in TifcZ with SageTj^X 

(The code in this section should work with any reasonable version of TifcZ, 
which means it should work with all but the most terribly out-of-date TeX 
installations — but to make sure we can accomodate everyone, the code here is 
commented out. You can almost certainly uncomment and run them. Make 
sure you do \usepackage-[tikz} in the preamble.) 



7 The sagecommandline environment 

When writing a T^rX document about Sage, you may want to show some ex- 
amples of commands and their output. But naturally, you are lazy and don't 
want to cut and paste the output into your document. "Why should I have to 
do anything? Why can't Sage and TjtX cooperate and do it for me?" you may 
cry. Well, they can cooperate: 



sage : 1+1 1 

2 2 

sage : is_prime (57) 3 

False 4 



Note that the output of the commands is not included in the source file, but 
are included in the typeset output. 

Because of the way the environment is implemented, not everything is ex- 
actly like using Sage in a terminal: the two commands below would produce 
some output, but don't here: 



sage: x = 2010; len (x . divisors () ) 5 

sage: print ' Hola , u mundo ! ' 6 
The difference lies in the Python distinction between statements and expres- 
sions. 

One nice thing is that you can set labels by using an @ sign: 

sage: 1 = matrix ([ [1 , 0 , 0] , [3/5 , 1 , 0] , [-2/5 , -2 , 1] ] ) 7 

sage: d = diagonal_matrix ( [15 , -1, 4]) 8 

sage: u = matrix ([ [1 , 0 , 1/3] , [0 , 1 , 2] , [0 , 0 , 1] ] ) # foo 9 

sage: l*d*u # this is a comment 10 
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[15 0 5] 11 
[9-1 1] 12 
[-6 2 6] 13 

And then refer to that label: it was on line ??, which is on page ??. Note 
that the other text after the hash mark on that line does not get typeset as 
a comment, and that you cannot have any space between the hash mark and 
the @. 

You can also typeset the output: 
sage : l*d*u 14 




sage: x = var('x') 18 
sage: ( 1 - cos (x) " 2 ) . tr ig_s impl if y ( ) 19 

sin (x) 2 

The Sage input and output is typeset using the listings package with the 
styles Sage Input and SageOutput, respectively. If you don't like the defaults 
you can change them. It is recommended to derive from Def aultSagelnput 
and Def aultSageOutput, for example. . . makes things overly colorful: 

sage : pi .n(100) 21 

3.1415926535897932384626433833 22 

Plotting things doesn't automatically pull in the plot, just the text repre- 
sentation of the plot (the equivalent of applying str() to it): 

sage: plot ( sin (x) , (x, 0, 2*pi)) 23 
Graphics object consisting of 1 graphics primitive 24 

You can include output, but it will be ignored. This is useful for doctesting, 
as all the sagecommandline environment things get put into the "_doctest . sage" 
file. 

sage: factor (x~2 + 2*x + 1) 25 
( x + 1 ) " 2 26 
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